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ABSTRACT 


We wish to receive a radar echo from the planets using the 


a 


facilities of the Arecibo Radio Observatory. * Calculations show that 


Venus, Mars, Mercury and Jupiter are likely radar targets for this 
facility. The planets have a predetermined motion with respect to the 


radar on the earth's surface, and to tune a narrow-band receiver to 
the frequency of the reflected echo, Doppler shifts must be predicted. 


An approximate equation for the relative velocity of a planet:is 

' 
PA = Pp 
geocentric co-ordinates of the planet. Plots of relative velocity, and 


V - 0.4413 cos 6 cos TN km/sec, where p, 6, @ are the 


Doppler shift are given, as well as plots of range and declination for 


the likely targets (planets) for the years 1960-1967. 


“NTRODUCTION 


T%., cadar under construction for the Arecibo Radio Observatory 
(ARO) may be used as an interplanetary radar with a maximum obser- 
vation time of 2 = hours per day on any planet in the declination in- 
terval -2° to 38°N. Because of the large antenna aperture (300-meter 
diameter), the sensitivity «f this radar is greater than that of Jodrell 
Bank by 24 db, other parameters being assumed equal. 

The spectrum of a radar echo from a planet is determined by 
the rate of rotation, the aspect of the pole, and by the type of surface. i 
Cohen discusses these effects and arrives at a parabolic spectrum 
(Figure 1) Doppler shifted from the transmitted frequency fy by an 


amount proportional to the radial velocity of the planet with respect 


This designation has been changed to Department of Defense Ionospheric 
Research Facility. 


@e 
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Figure41. Parabolic Spectrum of Echo from Rough Rotating Sphere. 


to the radar. The spectrum is centered about 


ie 


where Vis the radial velocity menticned. 

Using his approach and making certain assumptions, Cohen 
predicts the half-power bandwidth B to be of the order of one kilo- 
cycle for Mars and one-hundred cycles for Venus. 

It is evident from the calculations that follow, that the half- 
power bandwidth of the planetary echoes may be small compared to the 
Doppler shift fn fo 
hee echoes, then some prediction must be made in order to tune the 
receiver to the proper frequency band. The purpose of this paper is 


to predict the values of fy, for the planets from which radar echoes 


may be received. 


er es 


If a narrow-band device is to be used to receive 


RADAR EQUATION 
The zadar equation expresaing signal power to noise power is 


Ss P,A’e a 

N° FKTB4qhere io 
where P, is the transmitter power, 

go is the radar cross section of the target, 

is the effective antenna area, 
is the receiver noise figure, 
is Boltzmann's constant, 
is the ambient temperature of the receiver, 


is the bandwidth of the receiver, 


~ won my > 


is the wavelength, 
r is the range to the target, and 


a is the cable efficiency (transmit and receive). 


By putting the design parameters for the space radar into Equation (4) 
and making reasonable estimates for the others, we arrive at an ex- 
pression for signal-to-noise ratio in terms of planet size and distance 
from the earth. Using the parameters, | 
P= 2.5 x 10° watte 
‘A= 0.5 area of a 300-meter dish = 3.53 x 104 m@ 


c= Oana’, where a is the planet radius 


KS 70570) om 


a= 4 (3-db loss) , 
we arrive at the following equation: 


2 
r = 0.97 x 10°4 SS : (2) 


Tr 
where a and _r are in meters. 
: Z : ; : * 
The figure 0.4 wa™ asa radar cross section is an estimate. 
Since, for a smooth infinitely conducting sphere, o = ra’, the cross 


section may very well be nearer this latter figure, and the choice of 


0.4 na“ would therefore be conservative. 


OBSERVING TIME 


The radar is being designed to have a beam-swinging capability 
of 20° from the zenith in any direction. The beam width for \ = .70m 
is approximately 4/6° . A planet will traverse the stationary beam in 2/3 
minute. If the planet passes through the zenith (declination = 48°) the 
maximum observing time of 2 hours is available. Observing time per 
day will vary from a few minutes (a few beam widths) to the full 2 S 


hours, depending upon the planet's declination. Observation time versus 


declination for the ARO facility is plotted in Figure 2. 


RADAR TO PLANETS, SIGNAL-TO-NOISE RATIO 


Table I gives approximate planet radii and minimum ranges from 


Earth. It is of interest to calculate the signal-to-noise ratio on 


o = 0.074 ane from reference 2. 
moon 
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Table I. Planetary Constants 


Mean distance Rad:us Approximate Minimum flight 
Planet from sun (A. U.) (km) minimum distance time (minutes) 
from earth (A. U.) 


Venus 02723 6, 100 0.28 4,7 
Mars 1.524 3, 400 0.52 8.7 
Mercury 0. 387 ‘2, 400 0.61 10.2 
Jupiter 5, 203 70, 000 4,2 70 
Saturn 9.539 58,000 8.5 

Uranus 19.191 25,000 18 

Neptune 30. 070 27,000 29 

Pluto 3929 3,200 * 38 

Earth 1.000 6, 400 oes 


a per pulse basis of echoes from these planets using Equation (2). When 
Equation (2) is expressed in terms of radius a in kilometers and range 


r in astronomical units (A. U.), the result is 
(3) 


By replacing a in Equation (3) by the appropriate value, and by 
taking rz aa for each planet, we arrive at the signal-to-noise- 
ratios per pulse shown in Table II. 

It is evident from Table II that Venus is a target that requires 
no signal-to-noise improvement techniques and that the planets Mars, 


Jupiter, and Mercury represent likely targets with the aid of some 


Table Il. Signal-to-Noise Ratios as Function of Range and at Approximate 
Minimum Range. 


Planet ein Sirtr .) 
N N min 
Sc a ea ee 
4 
Venus 0.75/r 120 
4 
Mars 0.23/r 6 | 
4 
Mercury 0.12/r 0.8 
A 4 
Jupiter 100/r 0. 32 
4 
Saturn 67/r ROlZ 
4 -4 
Uranus (12. 5/r 1.2 x 10 
4 
Neptune 15/r Bix 1077 
4 -8 
Pluto 0.21/r 3 x 10 


signal-to-noise improvement techniques. Saturn, Uranus, Neptune 
and Pluto are not very likely targets because of their great distance 
from the earth. There is interest in these planets for passive meas - 
urements, however, and information on declination and range for 


Saturn, Uranus, and Neptune is included in this paper. 


RELATIVE VELOCITY 


To find an expression for the relative velocity of a planet with 


respect to a point on a rotating earth, we shall set the geocentric co- 
ordinate system up in such a way that the planet lies always in the y-z 
plane (see Figure 3). This gives the planet co-ordinates of (Pp : i $p) : 
The earth, then, rotates with respect to this co-ordinate system at the 


rate of approximately one revolution per day. 


Figure 3. Co-ordinate System. 


The co-ordinates of the fixed point on Earth A are py, Gy, oy 
where Pa and $, are constants; 6, is a function of time and is the 
hour angle modified by the change in right ascension of the planet. This 
means that the angular velocity of the earth's rotation must either be 


considered as a function of time or approximated as a constant. 


Of Venus, Mars, and Jupiter, the planet having the greatest day-to-day 
variation in right ascension is Venus. This day-to-day change in 196] is always 
less than six minutes. If the rotation period of the earth is taken as 23° 56™ : 
and d6,/dt is assumed constant in our co-ordinate system, we introduce a 
maximum error in the tangential velocity of the ARO of approximately 0.24 
per cent. 

We have obtained IBM cards from the Naval Observatory, which give 
day -to-day information on the planets through 1967. These cards contain ap- 
parent right ascension, declination, and true distance from the earth. This 
means that all the information on Pp and dp and their time derivatives is 
available. 

It may be shown (see Appendix) that the approximate equation for the 


relative velocity Vpa? accurate to approximately ane part in 104 18 


Va = Rp - 04413 cos 6 cos @,km/sec (9) 
where Pp is expressed in km/sec, and 6 is declination. In this equation the 


second term is of the order of one per cent of Pp maximum, for the planets of 


interest as radar targets. In the plots only the Pp term is considered. 
DOPPLER SHIFT 


The Doppler shift is given by 


Ae et eae Vea 


ie ie aot gee 
fe) fe) 
With fo = 430 Mc/s, the equation is 
af=-5x 10° Pp (A. U. /Day) cps 
The Doppier shift n.ay be obtained from the radial velocity plots by multi- 


plying the numbers appearing on the ordinate by minus 5 kilocycles, for 


example, see the radial plot for Mercury 1960. 


with respect to a point on the earth's surface V 


the following manner. 


APPENDIX 


To arrive at an expression for the relative velocity of the planet 


FP 


to rectangular co-ordinates (Figure 3),we find that 


and 


P(X, Y, Z) = (0, Pp sin dp, Pp cos $p) ; 


A(x, y, 2) = Pe cos 0, sing, » Pa sin@,, sin >, » Pa cos $y) ; 


A(x, y, Z) = (k, cos 6, , k, sin 6, , k5) 


A 


Converting the polar co-ordinates of P and A 


A? we may proc:ed in 


where k, = Pa sin by and k, = Pa cos >, - The distance between P and 


A is: 


PA 


a 
c= An Meerg ie aa?) 


2 P ; 2 2 
|(-#1 <2# 8) + (Op sings - ky sin 04] + ( Pp £08 bp -)*| 


rae shee Doe 2 A hee 2 
5 (sin 8, + cos ea) + k, + Pp (sin 8, + cos ep) 


2 
2 


~ 2Pp k, sing, sin@, - sin®, - 2pyk, co8 | 


2 PA ; : 
5 +k, + Pp - 2 Pb k, sindp sin6, ~ 2A, k, co #p| 
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The quantity Pp is the true distance to the planet, and the lowest 
value this parameter is likely to assume is approximately 0.3 A.U. for 
the closest approach of Venus. This amounts to approximately 4.5 x 10? km 
since 1 A.W, = 14.9504 x 107 em. 

When Venus is used as an example, typical maximum values for 
Pp , op and On are: 

bp +10 km/sec , 
op & 107! .rad/sec ; 
en = 1074 rad/sec 


Evaluating Equation (5) using the appropriate values given yields 


. 10(5 x 10' 46x10? ~2x 10°) +103+3x 10" a 
ae [107 +108 + 1016 4 1012 - 101] 2 


The second term in the numerator and all terms in the denominator 
except the ps may be neglected with an accuracy of about one part in 
104 The second two terms in the first term of the numerator are also 
effects of one part in 104 Simplifying Equation (5) by neglecting these 


terms, we get the appraximate equation 

t A A ° 
Vpa = Pp + dp (-k) cos 6p sin@, tk, sin bp) - oN (k) sin gdp cos 64) emt) 
but 4, = ; - § where 6 = declination; Then 


T : 
cos = cosi(— - 6) = sind 
p (5 - 6) 


sing, = sin (5 - 6) = cos 6 


and 


ip 


Putting the appropriate constants into Equation (7) and substituting yield 


Vv - 6'(2005cos 5 - 6052 sin6 sin@,) - 0.4443 cos 6 cos Oy, 3 


Pa ~ Pp 
(8) 


where Pp» bp» 5', and & may be read from day-to-day information 
on the cards; aN is restricted by pointing ability to 70° < IN < 110°. 
A > 14, and -0.34 < cos On < 0.34. This para- 


meter is a function of the time of day the measurement is made and must 


Therefore 0.94< sin6 


be left for local adjustment. The terms involving 6, are small correcting 
terms on the first-order term Pp : 

Using the maximum values of the parameters for Venus and placing 
them in Equation (8), we find that the second term is one part in 104 of 


Po max and the third term is one part in 10°, Therefore, to a good 


approximation, Equation (8) becomes 


Vv - 0.4443 cosicos6, km/sec . (9) 


ms i] 
PA a Pp A 
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